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Abstract 

In this paper we study the Euler-Poincare equations in M. N . We prove local 
existence of weak solutions in W 2 ' P (R N ), p > N, and local existence of unique 
£SJ . classical solutions in H k (R N ), k > N/2 + 3, as well as a blow-up criterion. 

For the zero dispersion equation(o = 0) we prove a finite time blow-up of the 
classical solution. We also prove that as the dispersion parameter vanishes, the 
weak solution converges to a solution of the zero dispersion equation with sharp 
rate as a — > 0, provided that the limiting solution belongs to C([0, T); H (M. )) 
with k > N/2 + 3. For the stationary weak solutions of the Euler-Poincare 
equations we prove a Liouville type theorem. Namely, for a > any weak 
solution u G H 1 (R N ) is u = 0; for a = any weak solution u G L 2 (R N ) is 
u = 0. 

Key words: finite time blow-up, zero dispersion limit, Liouville type theorem, Euler- 
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1 Introduction 

We consider the following Euler-Poincare equations in M. N : 

d t m + (u ■ V)m + (Vit) T m + (div u)rn = 0, 
(EP) { m = (1 - aA)u, 
u Q (x) = u , 



O 
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where u : R — > R is the velocity, m : R — >■ R represents the momentum, 
constant ^/a is a length scale parameter, (Vtt) T = the transpose of (Vu). The Euler- 
Poincare equations arise in diverse scientific applications and enjoy several remarkable 
properties both in the one-dimensional and multi-dimensional cases. 

The Euler-Poincare equations were first studied by Holm, Marsden, and Ratiu in 
1998 as a framework for modeling and analyzing fluid dynamics [13 [19], particularly 
for nonlinear shallow water waves, geophysical fluids and turbulence modeling. There 
are intensive researches on analogs viscous or inviscid, incompressible Lagrangian 
averaged models. We refer to [TJ H2J [26] for results on Navier-Stokes-a model in 
terms of existence and uniqueness, zero a limit to the Navier-Stokes equations, global 
attractor, etc. We refer to [2j [201 123] for results on analysis and simulation of vortex 
sheets with Birkhoff-Rott-a or Euler-a; approximation. 

For one-dimension, the Euler-Poincare equations coincide with the dispersion-less 
case of Camassa-Holm (CH) equation [I]: 

{CH) d t m + 2ud x m + d x u m = 0, m = (1 — ad xx )u . 

The solutions to (CH) are characterized by a discontinuity in the first derivative 
at their peaks and are thus referred to as peakon solutions. (CH) is completely 
integrable with a bi-Hamiltonian structure and their peakon solutions are true solitary 
waves that emerge from the initial data. Peakons exhibit a remarkable stability-their 
identity is preserved through nonlinear interactions, see, e.g. jU[22]. There are many 
comprehensive analysis on (CH) in the literature. We refer to a review paper [25] 
for a survey of recent results on well-poseness and existence of local and global weak 
solutions for (CH). The existence of a global weak solution and uniqueness was proven 
in [3l EJ [TUl El [29]. A class of the so called weak- weak solution was studied in [29] . 
The breakdown of the solution for (CH) was studied in [24]. 

The Euler-Poincare equations have many further interpretations beyond fluid ap- 
plications. For instance, in 2-D, it is exactly the same as the averaged template match- 
ing equation for computer vision (see, e.g., [TH [TTJ, [21] )• The Euler-Poincare equa- 
tions also has important applications in computational anatomy (see, e.g, |22, 30J). 
The Euler-Poincare equations can also be regarded as an evolutionary equation for a 
geodesic motion on a diffeomorphism group and it is associated with Euler-Poincare 
reduction via symmetry [Tj [TTJ, [221 E2 ED]- We refer to a recent book [22] for a 
comprehensive review on the subject. 

The organization of the paper is as follows. In Section 2, we give some preliminary 
discussions of the Euler-Poincare equations and we state a theorem on local existence 
of weak solution in W 2,P (W N ), p > N, and local existence of unique classical solutions 
in H k (R N ), k > N/2 + 3. 

In Section 3, we prove a theorem on a blow-up criterion, as well as, a theorem on 
finite time blow-up of the classical solution for the zero dispersion equation. For classic 
solutions with reflection symmetry, the divergence V • u satisfy a Riccati equation 
at the invariant point under the reflection transformation and hence there is a finite 
time blow up if the divergence is initially negative. 
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In Section 4, we prove that as the dispersion parameter a vanishes, the weak 
solution converges to a solution of the zero dispersion equation with a sharp rate 
as a — > 0, provided that the limiting solution belongs to C([0, T); H k {R N )) with 
k> N/2 + 3. 

Finally, for the stationary weak solutions of the Euler-Poincare equations we prove 
a Liouville type theorem in Section 5. For a > 0, we prove that any weak solution 
u G H 1 (M. N ) is u = 0. For a = 0, any weak solution u G L 2 (WL N ) is u = 0. This is a 
surprising result, as all the previous Liouville type results are for dissipative systems. 
This is the first Liouville type theorem for non-dissipative systems. 

We also give a proof of the local existence and uniqueness theorem in Appendix. 



2 Preliminaries and local existence 

In this section, we discuss some mathematical structures of (EP) and then we state 
a local existence theorem for the weak solution and the classic solution. We refer to 
[TTl 122] for more in-depth discussions on (EP). 
(EP) can be recast as 

d t m + V ■ (u(g>rn) + (Vu) T m = . (1) 

The last term above can be written in a conservative/tensor form 

N N N 

diUjirij = ^2 diUjUj — a ^2 diUjd 2 Uj 

3=1 i =1 hk=l 

I N N 

= -di\u\ 2 - a ^ d^diUjdkUj) + a ^ d k diUjd k Uj 

j,k=l j,k=l 
I N N 

= -di\u\ 2 - a ^ dj{diU k djU k ) +^^2 d i( d kUj) 2 

j,k=l j,k=l 



Set stress-tensor 



Tij = rriiUj + — \u\ 2 — adiu ■ dju + < ^^-\'Vu\ 2 . 



Then (EP) becomes 

N 

dtrrii + ^2 d 3 T a = • ( 2 ) 

3=1 

The first term in involves a second order derivative of u and it can be rewritten 

as 

N N 



rriiUj = UiUj 



+ a dkUidkUj - a^dk (ujd k Ui) . 



k=l k=l 
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The symmetric part of tensor T is given by 

T a = u®u + aVuVu 1 - aX7u T \7u + ^(\u\ 2 + a|Vu| 2 ) Id (3) 
and the remainder terms in T are given by 

N 

= -aJ2 d k( u Aui) ■ (4) 

fc=i 

Hence T = T a + T b . In view of this, the natural definition of the weak solution of 
(EP) would be: 

Definition 1 u G L°°(0,T; if* oc (R )) is a weak solution of (EP) with initial data 
u G Hl oc (R N ) if the following equation holds for all vector field <fi(x,t) such that 
</>(-, t) G Cq 00 ^) for all t G [0, T) and d>(x, ■) G Cft([0,T)) for all x G M. N 

[ [ (u- tf} t + aVu:V(j> t )dxdt+ [ («o • </>(•, 0) + aVu : V<£(-, 0)) dx 

Jo Jr n Jr n 

+ [ [ T a :Vdy(x,t)dxdt + a V" / / Ujdkmdjdkfadxdt = 0, (5) 

where T a is given by (T^j. 

(EP) also has a natural Hamiltonian structure. Set 

H = - I u ■ mdx. 

2 Jrn 

then = u and (EP) can be recast as 

d t m = -A^, (6) 
dm 

where A is an anti-symmetric operator defined by 

N N 

Au = dj{miUj) + djUjirij. 
i=i i=i 

Consequently, from (j2]) and ([6]), there are two conservation laws 

^- [ mdx = 0, — [ (\u\ 2 + a\Wu\ 2 ) dx = 0. 

For the one-dimensional case, (EP) coincides with the dispersion-less case of Camassa- 
Holm (CH) equation and there is an additional Hamiltonian structure and a Lax-pair 
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which leads to a complete integrability of (CH) [3]. We refer to [13] for a general 
discussion on bi-Hamiltonian system and complete integrability. 

When a = 0, the above Hamiltonian structure shows that (EP) is a symmetric 
hyperbolic system of conservation laws 



d t u + div(n <g> u) + -V\u\ 2 = 



(7) 



u(x, 0) = u 
which possess a global convex entropy function 

\d t \u\ 2 + divO|V) = 0. (8) 

We refer ([7]) as the zero dispersion equation. Indeed, we can recast it in a usual form 
of a symmetric hyperbolic system (we state it in 



u t + Au x + + Cu z = 

with 
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A is a symmetric matrix and has three eigenvalues: u, 2u + \u\, 2u—\u\, corresponding 
to one linearly degenerate field, and two genuinely nonlinear fields, respectively, when 
u ^ 0. 

We shall remark that although the high dimensional Burgers equation has a similar 
structure as ([7]), it does not possess a global convex entropy. In section 5, we will 
prove a Liouville type theorem for (J7J). This theorem does not hold true for the high 
dimensional Burgers equation. 

Now we introduce some notations and then we state a theorem on local existence 
of the weak solution and local existence and uniqueness of the classical solution. 

For s G R and p G [1, oo] we define the Bessel potential space L S,P (M. N ) as follows 

L S ' P (M. N ) = {/ G L P (R N ) | ||(1 - A)f f\\ LP := \\f\\ L s, P < oo}. 

For s G N U {0} it is well-known that L 8,P (M> N ) is equivalent to the standard Sobolev 
space W /s ' p (IR Ar )(see e.g. [2Z])- This, in turn, implies immediately that there exist 
Ci, C 2 such that 

Ci||u|| W /fe+2, P < ||m||£fc, P < C , 2|| , u||vi/fc+2,p (9) 
for all k G N U {0},j9 G (1, oo). As usual we denote H S (R N ) = W S > 2 {R N ). 

Theorem 1 (i) Assume a > and u G W 2 ' P (JH N ) with p > N. Then, there 
exists T = T(\\u Q \\w2,p) such that a weak solution to (EP) exists, and belongs to 
u G L°°(0, T; W 2 ' P (R N )) fl Lip(0, T; W 1 ' P (R N )). 
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(ii) Let a > and u G H k (R N ) with k > N/2 + 3. Then, there exists T = 
T(\\uo\\ H k) such that a classic solution to (EP) exists uniquely, and belongs to 
u G C([0, T); H k (R N )). 

(Hi) For a = 0, (EP) is a symmetric hyperbolic system of conservation laws with a 
convex entropy. Consequently, if Uq G H k (R N ) with k > N/2 + 1. Then, there 
exists T = T(||tto||.fffc) such that a classic solution to (EP) exists uniquely, and 
belongs to u G C([0, T); H k (R N )). 

The proof of symmetric hyperbolicity and existence of convex entropy in (iii) are 
given in fl7])-(JH]). The proof of existence of the unique classic solution for symmetric 
hyperbolic system is standard, see e.g [TH]. The proof of (i) and (ii) is also rather 
standard and will be given in the Appendix for completeness. 



3 Finite time blow up 

In this section, we first present a theorem on a blow-up criterion and then we prove 
a theorem on finite time blow up for the zero dispersion equation. 

We denote the deformation tensor for u by S = (Sij), where SV,- := \{diUj + djui). 
We recall the Besov space ^ which is defined as follows. Let {ip m }m& be the 
Littlewood-Paley partition of unity, where the Fourier transform i/> m (£) is supported 
on the annulus {£ G M. N \ 2™- 1 <\£\< 2 m }(see e.g. [28]). Then, 

/ e B* if and only if sup ||^ m * /||z~ : = ||/||b° < °°- 

mez 00,00 

The following is a well-known embedding result, 

L°°{R N )^BM0{R N )^B^ tO0 {R N ). (10) 

Theorem 2 For a > 0, we have the following finite time blow-up criterion of the 
local solution of (EP) in u G C([0,Q; H k (R N )), k > N/2 + 3. 

lim sup 1 1 it (t) l^* = oo if and only if / HS'WHbo dt = oo . (11) 
*-><» Jo 

Remark 1.1 Combining the embedding relation, W 1>N (M. N ) ^ BMO(R N ) ^ B^ j00 (R N ) 
with the inequality ||D 2 n|| L p < C||m|| iP for p G (l,oo)(see (IT5|) below), we have 

IISILo < C\\S\\ BM o < C\\DS\\ L N < C\\D 2 u\\ L N < C\\m\\ L N. 

Therefore we obtain the following criterion as an immediate corollary of the above 
theorem: for all p > N 

lim sup ||*n(i)||jy> = oo if and only if / \\m(t)\\ l n dt = oo. (12) 
*->t* Jo 
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Remark 1.2 In the one dimensional case of the Camassa-Holm equation (CH) the 
above criterion implies that finite time blow-up does not happen if J Q * Hw^r) H^i^t < 
oo for all t > 0. Thanks to the conservation law we have sup 0<r<< [[w^r)^ < 
1 1 ito || if 1 < 00 f° r ai l t > 0. Since we have embedding W /2 ' 1 (1R) ^ H 1 (M), and we do 
have finite time blow-up for (CH) [24], our criterion is sharp in this one dimensional 
case. 



Proof of Theorem 2 We only give a proof for the case a > 0. The proof for the 
case a = is similar and simpler hence will be omitted. 

Using estimates ( |33l 1341 1351 1361) for Ji,/ 2 ,^3 in the proof of Theorem 1 in the 
Appendix, one has 

— \\m(t)\\ H k < CdlVwIli^o + Umllioo + ||Vm||jr,oo)||m(t)||jr* 
< C(||m|| LP + \\Dm\\ LP + \\D 2 m\\ L v)\\m{t)\\ Hk . 

Hence, 

C [ {\\m(r)\\ LP + \\Dm(r)\\ LP + \\D 2 m(r)\\ LP }dr 
. Jo 

(13) 

for k > N/2 + 1 and p > N, where we used the Sobolev embedding. Consequently, 
blow up of ||ra(£)||#fc as t — > t* implies that at least one of ||m(t)||x, P , ||Z)m(t)||LP 
and ||Z^ 2 m(t)||iP blow up as t — > t* . In the following three steps, we show blow-up 
criterion for each of them are all given by (fTTj) . 

Step 1. We first recall the following logarithmic Sobolev inequality (see e.g. |28j), 

<c(i+imu )(io g (i + 11/11 (14) 



m {t)\\H k — \\ rn o\\H k ex P 



where s > 0,1 < p < 00 and sp > N. From the estimate in (1311) in the Appendix we 
obtain 

— ||m|| LP < C{1+\\S\\bo )log(l + ||S , || W ri J ,)||m|| LP (for p > N) 
at 

< C(l + II^II^^Jlog^ + llD^II^Hmll^ 

< C(l + ||5 , ||so o!0o )log(l + ||m|| LP )||m|| LP 

for p > N, where we used the boundedness on L P {R N ) of the pseudo-differential 
operator 

o- i:j (D) := didj{l - aA)- 1 = -RiRj A(l - aA)" 1 

with the Riesz transforms {Rj}^ =1 on R Ar (see Lemma 2.1, pp. 133[27J), which provides 
us with 

N 

\\D 2 u\\ LP = \Wij(D)m\\ LP < C\\m\\ LP (15) 
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for all p G (1, oo). By Gronwall's lemma we obtain 

log (1 + \\m(t) \\ LP ) < log(l + Hmollip) exp (c /(l + \\S(r)\\^ )dr) (16) 

V Jo / 

for p > N. This implies that 

limsup ||m(t)||x, P = oo if and only if / || -S'(i) || dt = oo . (17) 

t-+t* Jo 00 ' 00 

Step 2. Taking derivative of (EP) and L 2 (M. N ) inner product it with Dm\Dm\ p ~ 2 , 
we find that 

- — WmCtWr, = - / (div u)\Dm\ p dx - / (D« • V)m • ^mlDm^" 2 cfe 

P</t ' pJ R N J R N 

- D(Vu) T m- Dm\Dm\ p - 2 dx - (Vu) T Dm ■ Dm\Dm\ p - 2 dx 

Jr n Jr n 

— / D(div u)m ■ Dm\Dm\ p ~ 2 dx — / (divu) Dm ■ Dm\Dm\ p ~ 2 dx 

Jr n Jr n 

<|3 + -l/ \Du\\Dm\ p dx + 2 \D 2 u\ \m\ \Dm\ p - 1 dx 



< ^3 H — J ||Dn||ioo ||Dm||^ p + 2||D 2 n|| i 2 P ||m|| L 2 P llDmll^p 1 

< C||m||ip||Z)m||^p + CUmll^pH-Dmll^p 1 

for p > N, where we used the Sobolev embedding and ( |T5l) to estimate 

||Z)n||i,oo < C||-D 2, a||LP < (7||m||iP 
for p > N. Hence, for p > N we have 

— ||Dm(t)||ip < C||m|| LP ||Dm|| LP + C||m|| 2 2 p. 
By Gronwall's lemma, we have 

\\Dm(t)\\ LP < exp fc J" \\m(r)\\ LP dA (\\Dm \\ LP + C \\m{r)\\ 2 L2p dr^ (18) 

for p > N. ^From estimate (|16l) . one has 

f* \\m(s)\\ LP ds < tmax < s < ( \\m(s)\\ LP 
< tmax < s < t exp (log(l + ||m(s)|| LP )) 

< texp (log(l + ||m || L p)exp (c/ *(l + \\S(r)\\^ L Jdr^ . (19) 
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Similarly, 

J \\m(s)\\ L 2 P ds < texp ^log(l + ||ra |U 2 p) exp J (1 + \\S(t) \\^ J)drjj . 

(20) 

Combining (f!5 | PI5]) and (1201) . one obtains 

lim sup ||Dm(t) = oo if and only if / H^C^Ubo dt = oo . (21) 

Sfep 5. Similarly taking .D 2 of (EP) and L 2 (R N ) inner product it with D 2 m\D 2 m\ p ~ 2 , 
we find that 

~\\L^m(t)\\ p TP < 4 / |Dn||D 2 m| p rfx + 3 / |D 2 u| |Dm| |D 2 m| p " 1 da; 



+2 / \D 3 u\\m\\D 2 m\ p - 1 dx 



< 4||Dw|| L oo||£) 2 m||^ p + 3||D 2 n|| L 2 P ||L)m|| L 2 P ||D 2 m||^ 

+2||£) 3 w|| L 2 P ||m|| L 2 P ||D 2 m||^ 1 

< C||m|| i p||L' 2 m||^p + C||m|| L 2p||Dm|| L 2p||D 2 m||^ 1 
for p > N, where we used the estimate (TT5T) as follows 

||-D 3 u|| i9 = \\{D 2 (l-aA)- 1 }D(l-aA)u\\ Lq 

N 

< 2_. \\< J ij(D)Dm\\ Lq < C\\Dm\\ Lq , 

which holds for all q G (1, oo). Hence, 

^||L> 2 m(t)|| L p < C||m|| L p||D 2 m|| L p + C||m||j^ P ||I?m||i3i.. 
By Gronwall's lemma we have 

\\D 2 m(t)\\ LP < exp (c J \\m(T)\\ LP dr^ (\\D 2 m \\ LP + C \\m(r)\\ L 2 P \\Dm(T)\\ L 2 P dT 

for p > N. Similarly to the estimates in (fT9l) and fl20|) . the right hand side terms in 
the above inequality can all be controlled 

t : i +\\S(r)\\ B o)dr. 



Therefore, we have 

limsup \\D 2 m(t)\\ LP = oo if and only if / HS^H^o dt = oo . (22) 

Combination of ( TT3l 1T71 [2T1 [22}) gives the proof of the theorem. □ 
We now present a finite time blow-up result for a = 0. 
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Theorem 3 Let the initial data of the system ([?]), u G H k (R N ), k > N/2 + 2, has 
the reflection symmetry with respect to the origin, and satisfies divix (0) < 0. Then, 
there exists a finite time blow-up of the classical solution. 

Proof Taking divergence of (171), we find 

N N N 

d t {divu) + u-V{dwu) + 2 + ^{Au^uj + (divu) 2 + (didju^uj = 0, (23) 

where we used = \{diUj + djUi), and the fact 

N N N N N 



2j diUjdjUi + 2J diUjdiUj = 2 2J diUjSij = (SjUj + djUi)Sij = 2 2J S| 



^2 

i,j=l *)J = 1 *J=1 i , j = 1 i>i=l 



Now we consider the reflection transform: 

i? : x — )• x = —x, u(x,t) — >■ u(a;, t) = — it(— a;, i). 

Obviously the system (J7J) is invariant under this transform. The origin of the coor- 
dinate is the invariant point under the reflection transform. We consider the smooth 
initial data Mo G if (R ), A; > N/2 + 2, which has the reflection symmetry. Then, by 
the uniqueness of the local classical solution in H k (M, N ), and hence in C 2 (M, N ), the 
reflection symmetry is preserved as long as the classical solution persists. We consider 
the evolution of the solution at the origin of the coordinates. Then, u(0,t) = and 
D 2 u(0, t) = for all t G [0, T*), where T* is the maximal time of existence of the clas- 
sical solution in H k (R N ). If T* = 00, we will show that this leads to a contradiction. 
The system (123"]) at the origin is reduced to 



TV 

d t (divu) +2j2 S ij + (divw) 2 = 0, 

i,i=l 

which implies 

«9 t (div u) = -2 Si - (div n) 2 < -(div u) 2 . (24) 

Thus, 

divttn(O) 

dvu(O.i) < , . , 

v ' ' ~ l + divu (0)t 

which shows T* < | div u (o)\ ^ or divu (0) < 0. □ 

4 Zero a limit for weak solutions 

In this section, we show the following theorem on the zero dispersion limit a — > for 
the weak solutions. 
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Theorem 4 Let u a e L°°((0, T); if 1 (R JV )) 6e a weafc solution with initial data itg to 
(£PJ Q>0,flfi(iMe £°°((0, T); H k (R N )) fl Lzp((0, T); H 2 (M. N )), k > N/2 + 3, 
be the classic solution with initial data Uq to (EP) with a = 0, i.e., (7p. Then, we 
have 

sup < t < r {||u«(t) - u(t)\\ L 2 + vE||V(u a (t) - u(t))\M 

<C(a+ ||«« - w || L2 + v^HVK " ^o)IUO , (25) 

w/iere C = C(\\u\\ L c ^ >T;H k ( M.N^, \\u\\ Lip ( ,T;H 2 {RX))) is a constant. 

Proof We denote m := u — aAu. Then (it, m) satisfy (EP) with a truncation 
term as below 

d t m + div (it ® m) + (Vtt) T m = -a {Au ( + div(u ® Au) + {Vu) T Au} . (26) 

Subtracting f[2"61) from the first equation of (EP), and setting m := m a — m and 
u := u a — u. we find 

9 t m + div(tt <g> m) + div(tt <g> m) + div(it <g> to) + (Vit) T m + (Vu) T m + (V?x) T m 
= a {Au t + div(u <g> Au) + (Vit) T Ait} (27) 

Taking L 2 (M Ar ) inner product (J27J) with it, and integrating by part, and observing 

JR^ Jr-^ 

Jr^ Jr n 



we obtain that 
1 d 
2di 



JRJV JRJV Jr^ 

+a / [u ■ {Au t + div(u <g) An) + (Vtt) T Au}] rfx 



:= /i + J 2 + / 3 . 
We estimate 



N N 



Ii = — / djUj(uj — aAuj)uj dx — } / Uidi{uj — aAuj)uj dx 



Ji + j 2 



where 



AT JV 



Ji = — Vj / diUi\uj\ 2 dx + a y. / didkUi{dkUj)Ujdx 



i,j=l i,j,k=l ' 

N 

+a y. / diUi(dkUj)dkUjdx 
< C||^)|| c2 (||u||i 2 +a||Vu||| 2 ), 
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and 

N N 



i N r N r 

= — — y " / Uidi\uj\ 2 dx — a >, / dkUidi(dkUj)uj dx 

1,3=1 i,3,k=l 

-- V / Mj^l^Mjl 2 ^ 
i,j,k=l 

l N N 

= 77 / J \ diUi\uj\ 2 dx + a 2_, i didkUi{dkUj)Uj dx 

N r N f 

+a ^ / d k Ui(d k Uj)diUj dx + - V" / 9^19^1 2 

?,?,fe=l i,3,k=l 



< C||^)|| c2 (||«||| 2 + a||Vu||| 2 ). 

J2 = 2. / UidiUj(v,j — aAuj) dx 

i,3=l jRN 
N N 



y / UidiUjUj dx + a >^ / d k UidiUjd k Uj dx 

i,3=l jRN i,3,k=l jRN 
N 

+a / Uidid k Ujd k Uj dx 



i,3,i 

< CIKOHo-dlull^ + allVulli,). 
One can estimate I3 immediately as 

^3 < IMli 2 + a2C (\\ U \\Lip(0,T;H 2 (R N )) + \\ U \\l°° (0,T;H3 

Summarizing the above estimates, we obtain 

j t (\\u\\l 2 +a\\Vu\\l 2 )<C\\u^ 

which implies by Gronwall's lemma that 

llttlli, + q|| Vu||i 2 < C,{a 2 + l|w(0)||| 2 + a||Vu(0)||| 2 ) 

where constant C\ depended only on ||u|| up(o,T;H 2 (R N )) an d IMU°°(o,T;.h' 3 (r jv ))- This 
completes the proof of theorem. □ 
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5 Liouville type theorem for stationary solutions 

In this section, we prove a Liouville type theorem for stationary solutions. Recall 
that the stationary weak solution defined in Definition 1 reduces to 

Definition 2 u G H l (M. N ) is a stationary weak solution to (EP) on Mr, if the 
following holds 

N r N „ 

/ / {ujUj + aVuj ■ Vuj} djifj dx + a y / UjVui -Vdjcpidx 
+ Y [ !^l\u\ 2 -ad l u-d j u + ^\Vu\ 2 \d j i Pt dx = (28) 

~[j^ n I 2 2 J 

for i — 1,- ■ ■ ,N and for all c/> G C™(R N ) . 

Theorem 5 (i) Let u G -ff 1 (IR Ar ) be a stationary weak solution to (EP) with a > 0. 
Then, u = 0. 

(ii) Let u G L 2 (M. N ) be a stationary weak solution to (EP) with a = 0. Then, 
u = 0. 

Proof For a > 0, one can write ([28"]) in the following form, 

J2 T*d m dx+yr T$ jh d j d k <p i dx = Q, (29) 
where T£j is defined in and we recall here 



Tf, = UiU^ + aVui ■ Vuj + — \u\ 2 — adiU ■ djU H — IVkI 2 , 

J J J 2 2 



and 

Tf jk = aUjd k Ui. 



corresponding to Tf- in 



in rtlb 



Let us consider the radial cut-off function a G C^°(R N ) such that 





if \x\ 


1° 


if x 



a(|x|) 

and < o~(x) < 1 for 1 < |x| < 2. Then, for each R > 0, we define 



13 



Choosing ifi(x) = Xi(T R (x) in (1291) . we obtain 



N „ N N 

= T^a R (x)dx+Y / I? j x j a i a R (x)dx+J2 f^ ik d k a R (x)dx 

n „ N r 

+ E / E / T^x^a^dx 

= h + I 2 + h + h + h. (30) 

The hypothesis w G if 1 (IR iV ) implies that T G L 1 (M Ar ). Thus, we obtain 



\h\ < 4 / |T a ||x||V(x| da: < 2||V<t||l« / \T"\ dx -> 



i? 



as i? — ?• oo by the dominated convergence theorem. Similarly, J 3 , J 4 , J 5 — )■ as 
i? ->■ oo. 

Thus, passing i? — >■ oo in ( 1301 . we have 



iV 

= lim £ / na R {x)dx 
R -+°° ~[ Jr n 



which implies u = 0. This gives (i). 

For the case a = 0. All the terms involving a drop and (ii) holds true. This 
completes the proof the theorem □ 

We remark that the Liouville type results in Theorem 5 is rather surprising, as all 
the previous Liouville type results are for dissipative systems. For the Liouville type 
results for the dissipative systems, see, e.g. [5]. Theorem 5 is the first Liouville type 
theorem for non-dissipative systems. 
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Appendix: proof of Theorem 1 

Proof of Theorem 1 The proof of local existence part is standard, and below we derive 
the key local in time estimate of u(t) E L°°([0,T); W 2 ' P (R N )) n Lip(0,T; W 1 ' P (R N )). 

Id, 



— H m 



p dt 



\p 

I LP 



— - [u ■ V)|m| p dx + / djUimimj\m\ p 2 dx 
V Jr n Jr n 

''I.I — 1 

Jr n 



+ / (div u)\m\ p dx 



--\ [ TriS)\m\ p dx + V / |m| p - 2 cix 

/V JRJV JR N 

< C||S|| L oo||m||£ p < C||Vt*|| L =c||m||^ < C||m||£\ (31) 



and therefore 



^-||m|| LP < C||m|| LP . 
at 



We thus have the following estimate on L°°(0,T; W 2 ' P (R N )), 

where T = ^ Uo ^ - . In order to have estimate of u in Lip(0,T; W 1,P (R N )), we take 
L 2 (R N ) inner product (EP) with the test function V € W^'^QR^) for p > N. Then, 



/ dtm-i^dx = / m(u-V)if)dx— / m-Vu-t^dx 

Jr n Jr n Jr n 

\l 



< C||m||i,p||«||x,oo || VV'H ^ + CII^IIlpII Vu|| L oo || W|| ^ 



which provides us with the estimate, 

||<9t' u ||L°o(0,T;VK 1 ^(IR JV )) < C||9tH|U«'(0,T;W- 1 .P(R JV )) < C II m ll L°° (0,T);Lp 

Hence, for all < ii < t 2 < T we have 

ft 2 



||u(f 2 ) - < / ||^«(t)||wi,p<i* < C(t 2 ~ ti)\\m\\ 2 Loo{0T . LP{m)) . 

Jti 

Namely, 

||' w IUjp(0,T;W 1 >P(R JV )) ^ C|l m llz / ° o (0,T;LP(R JV ))- 

This gives (i). Next we prove local in time persistency of regularity for u(t) in 
H k (R N )) with k > N/2 + 3. Let (3 = (ft, • • • ,ftv) be the standard multi-index 
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notation with \/3\ — /3i + h (3n- Taking H k (R N ) inner product (EP) with m, we 

find 



We write 




:= Ji + J 2 , 



and using the standard commutator estimate, we deduce 

h < ^WD^u-V)™- (u-V)D p m\\ L 2\\D p m\\ L 2 

\P\<k 

< C(||Vt*||Loo||m|| H fc + ||u|| H fc||Vm||x,oo)||m||jj-* (34) 

< C(||it||jyjv/2+i+ e llmllj^-fc + Hm-ll^fc-a llTnrlljyjv/a+i+s) Umll^ffc (Ve > 0) 

< C\\mf Hk 

if k > N/2 + 1, where we used the fact u — (1 — «A) _1 m, and therefore H^H^ < 
||m||ij s -2 for all seR. 

J 2 = ir / ( u - V)\D a m\ 2 dx = ~ V / (div«)|D^m| 2 da; 

< C||Vtt||Loo||m||^- fc < C||m|| H iv/2-i+ e ||m||^* (Ve > 0) (35) 

< C||m||^ fc 

if > iV/2 — 1. The estimates of I2, ^3 are simpler, and we have 

h + h < ||(Vu) T m|| jH -fc||m|| // fc < (7(||Vit||i,oo||m||#* + ||u||#fc+i||ra||Lcx,)||ra||($6) 

< C(||m||^jv/2-i+ £ ||m||^fc + ||m||^fe-i||m|| i yjv/2+ e )||m|| JJ -fe 

< C\\m\\% k 

if k > N/2. Summarizing the above estimates, we obtain 

^HwiIIh* < C || rrx|| |^ fc 
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for k > N/2 + 1, which implies 



\u(t)\\ Hk < - C ^ H n k Vt G [0,T), where T ' 



1 — CUfiolli^fct C||iio||_ff fe 
where k > N/2 + 3. 

We now prove uniqueness of solution in this class. Let (u±, mi), (u 2 , m 2 ) two solution 
pairs corresponding to initial data (iti,o, wi^o), (1*2,0, m 2i0 ). We set u = ui — it 2 , and 
so on. Subtracting the equation for (u 2 , ra 2 ) from that of (1*1, mi), we find that 

d t m + div (ui <8> m) + div (m Cg> m 2 ) + (Vu^m + (Vn) T m 2 = 0. (37) 

Let p > N. Taking L 2 (M 7V ) product of (|37|) with m|ra| p ~ 2 , we obtain 

- — \\rn(t)\\ p T p = — (l — -) [ (div ui)\m\ p dx — [ (div u)m 2 ■ m\m\ p ~ 2 dx 
pat V PJ Jr n Jr n 

- [ («.V)m..m|m|^<b- / (Vt^m-mlmr* 
— / (/V u) T m 2 ■ m\m\ p ~ 2 dx 

< CdldivuiHioollmll^j, + ||Vw||L o||m 2 ||Lp||m,||^p 1 + ||u|| LP || Vm 2 || 

+ ||VMi||z,°o||ra||£ P + HVwIlL^ll^llLpHmll^p 1 ) 

< C(||tti|| H «i + ||it 2 || H 0ll m l 

for k > N/2 + 3. Hence, 



I v 

I LP 



||m(t)|| LP < ||mo||j>exp \ C J (\\ux{T)\\ H k + \\u 2 (T)\\ H k)dr 

This inequality implies the desired uniqueness of solutions in the class L l (0, T; H k (M. N )) 
with k > N/2 + 3. This gives (ii). The proof of (iii) was explained at the end of 
Section 2. This completes the proof of Theorem 1. □ 
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